INFINITE PERIODIC DISCRETE MINIMAL SURFACES 
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Abstract. A triangulated piecewise-linear minimal surface in Euclidean 3- 
space R 3 defined using a variational characterization is critical for area amongst 
all continuous piecewise-linear variations with compact support that preserve 
the simplicial structure. We explicitly construct examples of such surfaces that 
are embedded and are periodic in three independent directions of R 3 . 
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1. Introduction 

The goal of this article is to show existence of examples of discrete triply-periodic 
minimal surfaces that are modelled on smooth triply-periodic minimal surfaces. For 
each smooth minimal surface model we consider, we show how one can find a variety 
of corresponding discrete minimal surfaces. We restrict ourselves to discrete surfaces 
with a high degree of symmetry with respect to their density of vertices, and thus 
they have a highly discretized appearance. The advantage of such a restriction 
is two-fold: (1) we can give explicit mathematical proofs of minimality without 
relying on numerics, and (2) we can make changes in the symmetries that would 
not be allowed in the smooth case. In general one can also consider discrete minimal 
surfaces with finer triangulations, but considering only highly discretized examples 
is still sufficient to show existence of many discrete triply-periodic minimal surfaces 
corresponding to a single smooth one. For motivational purposes, we first briefly 
introduce smooth minimal surfaces. 



1.1. Smooth minimal surfaces. Soap films that do not contain bounded pockets 
of air are surfaces that minimize area with respect to their boundaries. Smooth 
compact minimal surfaces are mathematical models for soap films, because they 
are (by definition) surfaces that are critical for area with respect to all smooth 
variations that fix their boundaries. By computing the first derivative of area for 
a smooth variation of a general surface, one finds that minimal surfaces are also 
those whose two principal curvatures at each point are equal and opposite. Since 
the mean curvature at a point on a surface is the average of the principal curvatures, 
a minimal surface is then one for which the mean curvature is zero at every point. 

l 
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The simplest example of a minimal surface is the flat plane. Another example is 
the catenoid, which is a surface of revolution that can be parametrized by 

(1) { (cosh x cosy, cosh x sin y,x) gl 3 |iel,j/6 (0,2n] C R} , 

where R denotes the real numbers and R 3 denotes Euclidean 3-space. (Note that 
by restricting to soap films not containing bounded pockets of air, we have ruled 
out surfaces like the sphere, which is certainly a soap film but also has nonzero 
mean curvature.) 

As our goal is to study discrete minimal surfaces, with smooth minimal surfaces 
playing only a suggestive role, we will not go into more detail here about basic 
properties of the smooth case. (Some fine general introductions to smooth minimal 
surfaces are 0, 0, 0, HU, Bl, CE| and gOJ.) We will simply go directly to this 
definition: 

Definition 1.1. A smooth minimal surface in R 3 is a C°° immersion f : A4 — > 
R 3 of a 2- dimensional manifold M whose mean curvature is identically zero; or 
equivalently, the map f is critical for area with respect to all smooth variations 
compactly- supported in the interior of M.. 

A triply-periodic smooth surface is one that is periodic in three independent 
directions of R 3 , i.e. / and f + Vj have equal images for three independent constant 
vectors v\,V2, V3 in R 3 . There are a wide variety of smooth triply-periodic minimal 
surfaces, as can be seen by looking at papers of H. Karcher, K. Polthier, A. Schoen, 
W. Fischer and E. Koch g], HU, [T2|, [T3], EH], for example. We show a few 
examples here: the central surface in Figure was named the superman surface by 
W. Meeks ^1] (one special case of surfaces of this type is the Schwarz D surface) ; 
the second surface in the second row of Figure is a generalized Schwarz P surface 
(one special case of surfaces of this type is the original Schwarz P surface); the 
smooth Schwarz CLP surface is shown in the lower-right of Figure and the 
smooth triply-periodic Fischer-Koch surface is shown in the lower row of Figure lTTl 

1.2. Defining discrete minimal surfaces. Recently, finding discrete analogs of 
smooth objects has become an important theme in mathematics, appearing in a 
variety of places in analysis and geometry. So, minimal surface theory being a 
subject in geometry that relies heavily on analysis, it is natural to ask how to 
define a discrete analog of smooth minimal surfaces. But there is no single definitive 
answer; the definition one chooses would depend on which particular properties of 
smooth minimal surfaces one would wish to emulate in the discrete case. 

For example, a definition by A. Bobenko and U. Pinkall (2] uses discrete in- 
tegrable systems, in analogy to smooth integrable systems properties of smooth 
minimal surfaces (or, more accurately, of smooth surfaces with possibly-nonzero 
constant mean curvature). Their definition is good from the viewpoint of inte- 
grable systems, but does not yield discrete surfaces that are critical for area with 
respect to variations of their vertices. 

Here we will take a variational point of view, so we wish to consider area- 
critical discrete surfaces. We choose this definition: A discrete minimal surface 
in R 3 is a piecewise linear triangulated surface that is critical for area with respect 
to any compactly-supported boundary-fixing continuous piecewise-linear variation 
that preserves its simplicial structure, as defined in [21] and [2E| and SectionEJhere. 
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FIGURE 1 . Smooth and discrete superman surfaces. (All the computer 
graphics here were made using K. Polthier's Javaview software \'2'2\.) 

Although we will define discrete minimal surfaces as just above, we should re- 
mark that, even from within the variational point of view, this is not the unique 
choice of a definition. For example, a broader definition by K. Polthier [21] uses 
"non-conforming triangulations" , unlike the triangulations here which will all be 
conforming. Non-conforming discrete surfaces are those for which adjacent tri- 
angles are required only to intersect at midpoints of their boundary edges, not 
along entire edges. Then discrete minimality can again be defined variationally. 
This broader approach is useful for finding pairs of noncongruent isometric discrete 
minimal surfaces. Such pairs are called conjugate minimal surfaces in the case of 
smooth minimal surfaces. 

1.3. Constructing discrete minimal surfaces. Just as for smooth surfaces, a 
triply-periodic discrete surface is one that is periodic in three independent directions 
of M 3 . As noted at the beginning of this article, we will show that one can generally 
construct a variety of discrete examples modelled on a single smooth example. 
(Hence, in this sense, there exist an even wider variety of discrete triply-periodic 
minimal surfaces than there are smooth ones.) There are at least two possible ways 
to do this: 

• Method 1 vary the choice of simplicial structure of some compact portion 
of the surface, and 

• Method 2 vary the choice of rigid motions of M 3 that create the complete 
surface from some compact portion. 

To explain these two methods in more detail, imagine a compact portion M of 
a smooth triply-periodic minimal surface with piece-wise smooth boundary dM 
consisting of smooth curves 71,..., 7„. Suppose that each jj is either a straight line 
segment or a curve in a principal curvature direction of M that also lies in a plane of 
R 3 . (In the latter case we call jj a planar geodesic, since it is necessarily a geodesic 
of M.) A larger minimal surface M' is constructed from M by including images of 
M under 180° rotations about the lines containing linear and under reflections 
through the planes containing planar geodesic -fj . (The fact that the larger portion 
M 1 is still a smooth minimal surface can be shown using complex analysis, see JHI , 
|12j . |18| . |2()j . for example.) The larger portion M' again has a piece- wise smooth 
boundary dM' consisting of line segments and planar geodesies, so this procedure 
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FIGURE 2. Discrete and smooth Schwarz P surfaces. 



can be repeated again on M' . Repeating this procedure on ever-bigger pieces of 
the surface a countably infinite number of times, one builds the entire complete 
surface. M is often called a fundamental domain of the complete surface. 

For example, the minimal surface on the left-hand side of Figure ^ is a fun- 
damental domain M of a complete triply-periodic smooth minimal surface. The 
boundary dM contains the eight vertices 



for some given positive reals x,y,z > 0. Then dM is a polygonal loop consisting 
of eight line segments, from pj to Pj+i for j = 1, 2, 7, and finally from p$ to p\. 
One can construct the entire complete surface, as described above, using only 180° 
rotations about boundary line segments. A larger piece of this complete surface 
can be seen in the center of Figure ^ For general values of x,y,z, the resulting 
complete triply-periodic surface is a superman surface. When x = y this surface 
represents Schwarz' solution of Gergonne's problem (see and for more on 
this), and when x = y — \/2 ■ z this surface is the Schwarz D surface. 

To construct triply-periodic discrete minimal surfaces modelled on smooth su- 
perman surfaces, we first find a discrete version of the smooth fundamental domain 
M . One such example, found numerically using Java View software .22' , is on the 
right-hand side of Figure ^ However, there are many ways the simplicial structure 
of the discrete version can be chosen, and a number of them are shown in Figure 0] 
(The examples in Figure 0] all have coarser simplicial structures, and hence showing 
their minimality is manageable by direct computation without using a computer, 
as we will see in Section 0]) The corresponding complete triply-periodic discrete 
minimal surfaces are then constructed in the same way as in the smooth case, by 
180° rotations about boundary line segments. This variety of simplicial structures 
all based on one smooth superman surface is an example of using Method 1. 

To demonstrate Method 2, consider the left-hand minimal surface M in the 
second row of Figure Its boundary dM is two squares in parallel planes, where 



pi = (0,0,0), p 2 = (z,0,0), P3 
p 5 = (x,y,0), p 6 = (0,y,0), p 7 



(x,0,z), p 4 = (x,y,z), 
(0,y, z), p 8 = (0,0, z), 
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one square projects to the other by projection orthogonal to the planes. As in 
the previous example, 180° rotations about boundary line segments produces a 
complete triply-periodic smooth minimal surface, and a larger piece of this surface 
is shown just to the right of M in Figure [2] This example is a Schwarz P surface. 
The original Schwarz P surface is the special case that the distance between the 
two parallel planes containing dM is y2 times shorter than the length of each edge 
in dM, and is conjugate to the Schwarz D surface. (For later use, we mention that 
the lower-right surface M in Figure [21 is also a portion of a Schwarz P surface, now 
bounded by six planar geodesies, each lieing within one face of a rectanguloid with 
square base. The entire complete surface can be built from M solely by applying a 
set of reflections in equally-spaced planes parallel to the faces of the rectanguloid. 
The top one-fourth of M equals the bottom half of M. So one could choose either 
M or M as the fundamental piece for constructing a complete Schwarz P surface.) 

As with the superman surface, to create discrete analogs of the Schwarz P surface, 
one can apply Method 1 and choose amongst many different simplicial structures 
for this surface. Two such possibilities are the first two discrete minimal surfaces 
in the first row of Figure [2] The first one has 4 squares in parallel planes amongst 
its edge set, and the second one has 5 squares in parallel planes amongst its edge 
set. In fact, one can make examples with any number > 3 of such squares in its 
edge set, as we will see in Section \5\ In the third figure in the first row of Figure [2] 
a larger portion of a resulting complete discrete triply-periodic minimal surface is 
shown. 

But let us return to a demonstration of Method 2 on the Schwarz P surface. 
Let P be a plane perpendicular to the planes containing dM so that P also contains 
two disjoint boundary edges in dM . Consider the surface we would have if we first 
included a reflected image of M across P, and then created a complete surface by 
180° rotations about all resulting boundary line segments. A part of this surface 
is shown in the third figure of the second row of Figure [21 (this part is bounded by 
three boundary components, one square and two rectangles). In this part in Figure 
El there are three parallel dotted lines, along which the surface is not even a C 1 
immersion, hence the mean curvature is not defined there, and we cannot call this a 
minimal surface. So this construction is forbidden in the smooth case. However, for 
the discrete analogs shown in Figure [21 such a construction actually does produce 
a discrete complete triply-periodic minimal surface. This construction works in the 
discrete case because the notion of " minimality" is defined only at the vertices of 
the discrete surface, as we will see in Section [21 Choosing the simplicial structure 
so that there are no vertices in the interiors of the edges comprising dM is what 
makes this construction possible. 

Thus we have varied the rigid motions of R 3 that were used to create the complete 
smooth surface from M, in a way that cannot be allowed for smooth surfaces, to 
make a different kind of discrete minimal surface modelled on M. Part of this 
surface is shown in the upper-right of Figure [21 (this part also has three boundary 
components, one square and two rectangles). One can easily imagine including more 
reflections (not just across a single plane P), to make infinitely many different kinds 
of discrete minimal surfaces modelled on M. This is Method 2. 

In Section we state our main result. In Sections 01 and El as applications of 
Method 1, we will see a variety of different simplicial structures making discrete 
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FIGURE 3. At each vertex p the gradient of discrete area is the sum of 
the 90°-rotated edge vectors J(r — q), as in Equation 

analogs of the smooth superman and Schwarz P surfaces. Section contains exam- 
ples modelled on other smooth minimal surfaces. In the last example of Section [5] 
and the first example of Section [S] there are further applications of Method 2. 

2. Discrete Minimal Surfaces 

In the smooth case, a compact minimal surface is area-critical for any variation 
that fixes the boundary. We wish to define discrete minimal surfaces so that they 
have the same variational property for the same types of variations. We begin 
by defining discrete surfaces and their variations, and we first give an informal 
definition: 

Definition 2.1. (Informal Definition) A discrete surface in R 3 is a C° mapping 
f : A4 — > R 3 of a 2 -dimensional manifold A4 so that each face of some triangulation 
of Ai is mapped to a triangle in R 3 . The surface f(M) is embedded if f is an 
injection. 

We define embedded in the discrete case without any conditions about nondc- 
generacy of / (nondegeneracy is meaningless here, as / is only C ). However, we 
still use this word, to maintain the analogy to embeddedness of smooth surfaces. 

We stated the above definition for the intuition it provides, but we will require 
a more involved definition, which we now give: 

Definition 2.2. (Formal Definition) A discrete surface in R 3 is a triangular mesh 
which has the topology of an abstract 2-dimensional locally-finite simplicial surface 
K combined with a geometric C° realization T in R 3 that is piecewise-linear on 
each simplex. (Because K is a simplicial "surface", each l-dimensional simplex in 
K lies in the boundary of exactly one or two 2-dimensional simplices of K .) The 
geometric realization T is determined by a set of vertices V = {pi,p2,...} C R 3 
corresponding to the ^-dimensional simplices of K , and V could be either finite or 
countably infinite. The simplicial surface K represents the connectivity of T . The 
0, 1, and 2 dimensional simplices of K represent the vertices, edges, and triangles 
ofT. 

Let T = (p,q,r) denote an oriented triangle of T with vertices p,q,r £ V. Let 
pq denote an edge of T with endpoints p,q<E V. 

For p £ V, let star(p) denote the triangles ofT that contain p as a vertex. 

The boundary dT ofT is the union of those edges bounding only a single triangle 
ofT. The interior vertices (respectively, boundary vertices J ofT are those that are 
not contained (respectively, are contained) in dT. 
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We say that T is complete if dT is empty and if T is complete with respect to 
the distance function induced by its realization in R 3 . 

Definition 2.3. Let V = {pi,p2, ■■■} be the set of vertices of a discrete surface T. 
A variation T(t) of T is defined as a C°° variation of the vertices pi 

Pi(t) : [0, e) — > R 3 so that Pi(0) = p% Vi = 1, m. 

The straightness of the edges and the flatness of the triangles are preserved as the 
vertices Pi(t) move with respect to t. 

When T is compact, we say that T(t) fixes the boundary dT if Pi(t) is constant 
in t for all pi € dT. When T is complete, we say that T{t) is compactly supported 
if Pi(t) is constant in t for all but a finite number of vertices pi. 



The area of a discrete surface is 

E 

TeT 

where area T denotes the Euclidean area of the triangle T as a subset of ' 



areaT := areaT , 



Lemma 2.1. Let T(t) be a variation of a discrete surface T. At each vertex p of 
T , the gradient of area is 

(2) V p areaT = i £ J{r - q) , 

T—(p,q,r) Gstar(p) 

where J is 90° rotation in the plane of each oriented triangle T. The first derivative 
of the area is then given by the chain rule 

(3) -^areaT(i) 



/ d{p{t)) 

=0 „ e y \ 



dt 



, V p areaT 

t=0 



Proof. Let Pi{t) be the corresponding variation of each vertex in the vertex set V{t) 
of the variation T(t). Then 



area(T(t)) = i ]T £ ||(r(i) - p(t)) x (g(i) - p{t))\\ 

p(t)€V(t) \(p(t), g (t),r(t))estar(p(t)) 

and a computation implies 

|area(T(i)) = i £ /», £ \\r(t) - q (t) M t)) , 

where is the unit conormal in the plane of the triangle (p(t),q(t),r(t)) along 
the edge r(t) — q(t), oriented in the same direction as J(r(t) — q(t)). Restricting to 
t = proves the lemma. □ 

As defined in Section^] a smooth immersion / : M. — > R 3 of a 2-dimensional com- 
plete manifold without boundary is minimal if / is area-critical for all compactly- 
supported smooth variations. In the case that M is compact with boundary, then 
/ is minimal if it is area-critical for all smooth variations preserving f(dM.). 

We wish to define discrete minimal surfaces T so that they have the analogous 
properties, for variations as in Definition 12.31 So when T is compact, we consider 
variations T(i) of T that fix dT; and when T is complete, we consider variations 
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T(i) of T that are compactly supported. By Lemma l2~Tl the condition that makes T 
area-critical for any variation of these types is expressed in the following definition. 

Definition 2.4. A discrete surface is minimal if 

(4) V p areaT = 

for all interior vertices p. 

Remark 2.1. If T is a discrete minimal surface that contains a discrete subsurface 
T' lieing in a plane P, it follows from Equations ||2J and Q that the discrete 
minimality of T is independent of the choice of triangulation of the trace of T 
within P. Thus whenever such a planar part T' occurs in the following examples, 
we will be free to triangulate T' any way we please, within its trace in P. 

3. Results 

For the purpose of stating our main theorem, we give the following two defini- 
tions: 

Definition 3.1. A discrete triply-periodic minimal surface T has common topology 
and symmetry as a smooth triply-periodic minimal immersion f -.M^M 3 if there 
exists a homeomorphism 

<j>:f{M)^T 

such that the following statement holds: R s : R 3 — > R 3 is a rigid motion preserving 
f{M.) if and only if there exists a rigid motion R c i : R 3 — * R 3 preserving T so that 

R d ocj) = <t>° R s \f( M ) , 

and furthermore R s is a reflection (resp. translation, rotation, screw motion) if 
and only if Rd is a reflection (resp. translation, rotation, screw motion). 

Definition 3.2. We say that a subsurface T' of a complete discrete triply-periodic 
minimal surface T is a fundamental domain if T' can be extended to all of T by a 
discrete group of rigid motions {Rd,a}a£A generated by 

(1) reflections across planes containing boundary edges and 

(2) 180° degree rotations about boundary edges 

so that each Rd, a is a symmetry of the full surface T . 

Remark 3.1. In the above definition of a fundamental domain, we do not allow 
rigid motions that do not fix any edges of T (thus any fundamental domain of 
the example in Subsection \4-l\ must contain at least 6 triangles, even in the most 
symmetric case x = 1). Also, we do not allow rigid motions that are not symmetries 
of the full surface T ( thus any fundamental domain of the example in Subsection 
\5.<H must contain at least 32 triangles). 

We now state our results about embedded triply-periodic discrete minimal sur- 
faces, which involve comparisons to the following smooth minimal surfaces: the 
superman surfaces (Figure [TJ, the Schwarz P surfaces (Figure HJ, the Schwarz H 
surfaces, the Schwarz CLP surfaces (Figure [TJl, A. Schoen's I-Wp and F-Rd and 
H-T surfaces, and the triply-periodic Fischer-Koch surfaces fFigure lTT|l . More com- 
plete information about these smooth surfaces can be found in 0], 0, ^Hli 

dSj, IH, ESI, EH, CHI, E3 and |2H1 
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Theorem 3.1. The following discrete embedded triply-periodic minimal surfaces 
exist: 

(1) those with common topology and symmetry as smooth superman surfaces 
whose fundamental domains contain 4, 5, 6 or 8 triangles; 

(2) those with common topology and symmetry as smooth Schwarz P surfaces 
whose fundamental domains contain 1, 2, 6 or 32 triangles, and also a 
different class of discrete surfaces with common topology and symmetry as 
smooth Schwarz P surfaces whose fundamental domains contain 2n trian- 
gles for any positive integer n; 

(3) those with common topology and symmetry as smooth Schwarz H surfaces 
whose fundamental domains contain 2n triangles for any positive integer n; 

(4) those with common topology and symmetry as smooth Schwarz CLP surfaces 
whose fundamental domains contain 6 triangles; 

(5) one with common topology and symmetry as A. Schoen's smooth I-Wp sur- 
face whose fundamental domain contains 5 triangles; 

(6) one with common topology and symmetry as A. Schoen's smooth F-Rd sur- 
face whose fundamental domain contains 3 triangles; 

(7) those with common topology and symmetry as A. Schoen's smooth H-T 
surfaces whose fundamental domains contain 6 triangles; 

(8) those with common topology and symmetry as the smooth triply-periodic 
surfaces of Fischer-Koch whose fundamental domains contain 8 triangles. 

To prove this theorem, we need only collect the examples proven to exist in the 
remainder of this paper, as follows: 

Proof. Embedded discrete superman surfaces whose fundamental domains contain 

4 (resp. 5, 6, 8) triangles are given in Subsection l4.3l fresp. 14.41 14.21 14.1JI . A second 
type of embedded discrete superman surfaces whose fundamental domains contain 6 
triangles are given in the second to the last paragraph of Subsection lG.ll bv Method 
2. 

Embedded discrete Schwarz P surfaces whose fundamental domains contain 1 
(resp. 2, 6, 32) triangles are given in the first example of Subsection 15 .11 (resp. the 
second example of Subsection 15. II Subsection 15. 21 the last paragraph of Subsection 
16. II by Method 2). Other classes of embedded discrete Schwarz P surfaces are given 
by choosing k — 4 and zq — and jo = n in Subsection 15.31 for any positive integer 
n, and here we can choose the fundamental domains to be the 2n triangles between 
two adjacent meridans and below the plane {{x, y, 0) | x, y S M}. 

Embedded discrete Schwarz H surfaces are given by choosing k = 3 and zq = 
and jo — n in Subsection l5 . 3l for any positive integer n, and here again we can choose 
the fundamental domains to be the 2n triangles between two adjacent meridans and 
below the plane {(x, y,Q) | x,y € R}. 

Embedded discrete Schwarz CLP surfaces whose fundamental domains contain 
6 triangles are given in Subsection 16. II 

Embedded discrete I-Wp and F-Rd surfaces whose fundamental domains contain 

5 and 3 triangles, respectively, are given in Subsection 16. 21 

Embedded discrete H-T surfaces whose fundamental domains contain 6 triangles 
are given in Subsection 16. 31 

Embedded discrete triply-periodic Fischer-Koch surfaces whose fundamental do- 
mains contain 8 triangles are given in Subsection 16.41 □ 
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4. Discrete versions of the superman surface 

In Sections 01 El and El we construct discrete triply-periodic minimal surfaces. 
All of the surfaces we construct are embedded. 

To construct examples, we always start with a compact discrete fundamental 
piece T, with given simplicial structure and boundary constraints. The complete 
triply-periodic discrete surface is then formed by including images of T under a 
discrete group of rigid motions of R 3 . This group of rigid motions is generated by 
a finite number of 180° rotations about lines and/or reflections across planes, and 
for each edge pq in dT this group contains either 

• the 180° rotation about the line containing pq, or 

• a reflection across a plane containing pq. 

To ensure that the resulting complete discrete triply-periodic surface is minimal, 
Section gave us the following two approaches: 

(1) Use symmetries of T and of the resulting complete discrete surface to show 
that Equation (J2J holds at the vertices. 

(2) Locate the vertices of T so that T is area-critical with respect to its bound- 
ary constraints. 

In the following examples, either approach produces the same conditions for mini- 
mality. 

As noted in the introduction, we wish to show examples here for which explicit 
mathematical proofs of minimality are still manageable (without the aid of a com- 
puter) . So we are limited to examples with a high degree of symmetry with respect 
to their density of vertices, and thus with a highly discretized appearance. Discrete 
minimal surfaces that appear more like approximations of smooth minimal surfaces 
usually can only be found numerically. Numerical examples, with finer simplicial 
structures, of discrete versions of the superman, Schwarz P, F-Rd, I-Wp and H-T 
surfaces are shown in |24|. 

4.1. First example. The fundamental piece T here has eight boundary vertices 

Pl = (1,0,0), P2 = (1,1,0), pa = (1,1,3!), p 4 = (0,1,2:), 

p 5 = (0,1,0), p 6 = (0,0,0), P7 = (0,0,x), pb = (1,0,3;), 
for any given fixed x > 0, and has one interior vertex 

Pv = (h 2' f) • 
There are eight triangles in T, which are 

(pj,Pj+i,P9), 3 =1, 7, (p 8 ,pi,p g ). 

The complete triply-periodic surface is generated by including the image of T under 
180° rotations about each edge of dT, and then continuing to include the images 
under 180° rotations about each resulting boundary edge until the surface is com- 
plete. It is evident from the symmetries of this surface that Equation Q holds at 
every vertex. The fundamental piece T and a larger part of the resulting complete 
surface are shown on the left-hand side and center of the first row of Figure 0] for 
x = 1. The case for some given x € (0, 1) is shown on the right-hand side of the 
first row of Figure 01 
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FIGURE 4. Four different discrete versions of the superman surface. 

4.2. Second example. The fundamental piece T here has six boundary vertices 
pi = (0,0,0), p2 = (x,0,Q), P3 = (x,y,0), 
p 4 = (x,y,l), p 5 = (0,y,l), p 6 = (0,0,1), 

for any given fixed x, y > 0, and has one interior vertex 

fi \ 2 ) 2 ' 2' ' 
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There are six triangles in T, which are 

{Pj,P] + l,P7), 3 — 1 , ■ ■ ■ , 5 , (p 6 ,P!,p 7 ). 

The complete triply-periodic surface is generated by 180° rotations about bound- 
ary edges, just as in the previous example. In this example as well, it is evident 
from the symmetries of this surface that Equation holds at every vertex. Two 
fundamental pieces T of different sizes and larger parts of the resulting complete 
surfaces are shown in the second row of Figure 0] (x = y = 1 in the first case, and 
x < 1 < y in the second case). 

In the case that x = y = 1, this fundamental piece T has the same boundary as 
a fundamental piece of the smooth Schwarz D surface. Furthermore, for general x 
and y, this surface can be viewed as a discrete analog of the superman surface as 
follows: Consider the eight-straight-edged polygonal curve from the point (0, 0, — |) 
to the point (x, —y, — |) and then to (x, —y, ^) and then to (2x, 0, |) and then to 
(2cc,0, — |) and then to (x, y, — i) and then to {x,y, i) and then to (0,0, i) and 
then back to (0, 0, —5). This polygonal curve is contained in this discrete surface 
(although not in its edge set) and is also the boundary of a smooth superman 
surface. 

4.3. Third example. The fundamental piece T here has four boundary vertices 

Pl = (0,0,0), P2 = (1,1,0), Pa = (l,l,z), P4 = (1,0,Z), 
for any given fixed z > 0, and has one interior vertex 

p 5 = (a, b, c) . 
There are four triangles in T, which are 

(j?j,Pj+i>P5), 3 = l,-..,3 , (P4,P1,P5). 

Reflecting T across the plane {(x%, 0, £3) € R 3 | xi, X3 el} and attaching its image 
to T, one has a larger discrete surface containing eight triangles and six boundary 
edges. One can extend this larger discrete surface to a complete triply-periodic 
surface by 180° rotations about boundary edges, just as in the previous examples. 
In this example, Equation Q holds at each vertex j?i,...,p4 in the resulting complete 
surface. However, getting this to hold at p$ requires proper choices of a and b. 

For simplicity, we restrict to the case 2 = 1. Then, by symmetry, we may assume 
b — c. A computation shows that Equation J3J holding at p$ is equivalent to 

(5) (1 - a)vV - 2ab + 3b 2 = (a - b)y/(l - a) 2 + (1 - bf , 



(6) (1 - b)^a 2 - 2ab + 36 2 = (36 - a)y/(l - a) 2 + (1 - b) 2 . 

The solution to this is 

(7 b = — , a — . 

y ' 2 ' 2 

So when z = 1 and b = c and Equation Q holds, the area gradient is zero at each 
vertex pj for j = 1, 2, 9 in the extended complete triply-periodic discrete surface, 
and then symmetries of the surface imply the entire complete surface is minimal. 

Since the above minimality condition JSJ-lO is a system of two equations in 
two variables a and 6, we say the minimality condition here (when z = 1) is two- 
dimensional. 
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FIGURE 5. Three different discrete versions of the Schwarz P surface. 



The fundamental piece T with z = 1 is shown on the left-hand side of the third 
row of Figure 01 and a larger part of the resulting complete surface is shown just 
to the right of this. 

4.4. Fourth example. The fundamental piece T here has five boundary vertices 

pi = (0,0,0), pa = (1,1,0), p 3 = (l,M), P4 = (0,M), p 5 = (0,0,*) 
for any given fixed z > 0, and has one interior vertex 

p 6 = (a, 1 - a, b) . 
There are five triangles in T, which are 

(PitPj+uPe), 3 = 1 , -A, (P5,Pi,Pe)- 

The complete triply-periodic surface is generated by 180° rotations about boundary 
edges. In this example, Equation Q holds at each vertex pi,...,p 5 in the resulting 
complete surface, and making it hold also at pe requires proper choices of a and b. 

Like in the previous example, we can find a pair of explicit equations, in the 
variables a and b, that represent the minimality condition. These equations are 
similar to those of the previous example, and are slightly more complicated. One 
can then show the existence of a and b solving this minimality condition. 

Two fundamental pieces T of different sizes (z = 1 in the first case, and z < 1 in 
the second case) are shown on the left and right-hand sides of the bottom row of 
Figure 0] A larger part of the resulting complete surface in the case z — 1 is shown 
in the bottom-middle of Figure 0] 
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5. Discrete versions of the Schwarz P surface 

5.1. First two examples. Consider the vertices 

Pl = (3,0,6), p 2 = (6,0,3), p 3 = (6,3,0), 

p 4 = (3,6,0), p 5 = (0,6,3), p 6 = (0,3,6), p 7 = (3,3,3), 
and let 7i be the planar fundamental domain with the six triangles 

(Pj,Pj+i,P7), 3 — 1 , ■ ■ ■ , 5 , (p 6 ,Pi,p 7 )- 
Also, consider the vertices 

px = (3,0,6), p 2 = (4,0,4), p 3 = (6,0,3), p 4 = (6,2,2), 
p 5 = (6,3,0), p 6 = (4,4,0), p7 = (3,6,0), p 8 = (2,6,2), 
p 9 = (0,6,3), pio = (0,4,4), p u = (0,3,6), p 12 = (2,2,6), p 13 = (3,3,3), 
and let 72 be the fundamental domain with the twelve triangles 

(Pj,Pj+l,Pl3), j = 1, 11 , (Pl2,Pl,Pl3)- 

We can extend Tj (for either j = 1, 2) to a complete discrete surface by includ- 
ing the images of Tj under the reflections across the planes {(x,y,6k) \ x, y £ R}, 
{(x, 6fc, z) | x, z G K} and {(6fc, y, z) | y, z E R} for all integers k. Furthermore, 
Equation |@} holds at all vertices of the extended surface, so it is minimal. (See the 
first two columns of Figure |S] ) 

The surface produced by T\ (resp. 7^) is a simpler (resp. more complicated) 
version of a discrete Schwarz surface. Note that they are analogous to the bottom- 
right picture in Figure El (The second example T 2 was also shown in |25|.) 

5.2. Third example. Consider the ten vertices 

Pj = (a, (-l)'o, 1), p j+2 = (a, 1, (~iy +1 a), p j+i = (a, (-l) j+1 a, -1), 

p j+e = (1, (-iya, a), p j+8 = (1, (-iy +1 a, -a), j - 1, 2, 

and let T be the discrete surface with the eight triangles 

(Pl,P2,P7), (P2,PS,P7), (P2,P3,P»), (P3,P4,Ps), 
(P4,p 9 ,ps), (p4,p 5 ,p 9 ), (p 5 ,p 6 ,p 9 ), ( P 6,Pl0,P9)- 

Then let T be the discrete surface, with 24 vertices and 32 faces, that is made by 
including the four images of T under the rotations about the axis {(0, 0, r) \ r G M} 
of angles 0°, 90°, 180° and 270°. This T is shown in the upper-right of Figure 

One can then generate a complete triply-periodic surface by including the images 
of T under reflections across the planes {{x, y,k)\x,y € M}, {(x, k,z)\x,z 6 M} 
and {(k, y,z)\y 7 z £ R} for all integers k. The result of applying one such reflection 
is shown in the lower picture of the right-most column of Figure 

The condition for this discrete triply-periodic surface to be minimal is that 

_ 3V2-V3 

i.e. for this value of a, Equation (0} holds at every vertex of the surface. 
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5.3. Examples based on discrete minimal catenoids. Here we give two closely- 
related types of examples based on discrete minimal catenoids. One type is a dis- 
crete analog of the Schwarz P surface. The other type is actually an analog of the 
smooth Schwarz H surface, not the Schwarz P surface. To construct these exam- 
ples, we will use discrete analogs of the catenoid |25| . which are described in terms 
of the hyperbolic cosine function, just as the smooth catenoid was in Equation IjTJl. 

The vertices of a discrete minimal catenoid lie on congruent planar polygonal 
meridians, and the meridians are contained in planes that meet along a single line 
(the axis) at equal angles. Every meridian is the image of every other meridian by 
some rotation about the axis. By drawing edges between corresponding vertices of 
adjacent meridians (i.e. so that these edges are perpendicular to the axis), we have 
a piecewise linear continuous surface tessellated by planar isosceles trapezoids. We 
can triangulate each trapezoid any way we please without affecting minimality, as 
noted in Remark 12.11 so we shall triangulate each trapezoid by drawing a single 
diagonal edge across it. 

Two examples of discrete catenoids are shown in the first two pictures in the 
upper row of Figure El Both of these pictures have adjacent meridians in planes 
meeting at 90° angles. The first (resp. second) one has four (resp. five) vertices 
in each meridian. Another example is shown in the left-most picture of Figure El 
where the adjacent meridians lie in planes meeting at 120° angles, and there are 
four vertices in each meridian. 

To explicitly describe discrete catenoids, we need only specify: 

(1) The axis I: let us fix i = {(0, 0,z)\z€ K}. 

(2) The angle between planes of adjacent meridians: let us fix 8 — ^ for 
some integer k > 3. 

(3) The locations of the vertices along one meridian. 

We can place one meridian in the plane {(x, 0, z) \ x, z G K}, and locating its vertices 
at the following points will ensure minimality of the surface (see |25|1: 

pj = (rcosh ( -a(z + jS)) ,0,z +jS) 



with j = j , j + 1, ... , for some integers jo and j\ {jo < ji), and with 



r i lS 2 



= — arccosh 1 



5 V r 1 + cos < 



where r > and 6 > and zq 6 K are constant. The edges along this meridian are 
PjPj+i for j between jo and j\ — 1. 

For our application, we shall restrict to either k = 4, as in Figure[3 or to k = 3, 
as in Figure We shall further assume that cither 

• zq = and jo = —j\ < 0, or 

• zo = f and j = -ji - 1 < -1. 

Either of these conditions will produce a discrete minimal surface T whose trace 
has dihedral symmetry. One can then extend T by 180° rotation about boundary 
lines to a complete embedded discrete surface in R 3 . To conclude minimality of this 
complete surface, it remains only to check that Equation (0J holds at any vertex 
contained in any edge about which a 180° rotation was made, and this is clear from 
the symmetry of the surface. 
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FIGURE 6. Discrete version of the Schwarz H surface. 

The case when k — 4 and zq = and jo = —j\ = — 2 is shown in the second 
picture of the first row of Figure |3 and a larger portion of the resulting complete 
minimal surface is shown in the picture just to the right of it. The case when k = 3 
and Zq = | and jo = —j\ — 1 = —2 is shown in the left-most picture of Figure 
El and a larger portion of the resulting complete minimal surface is shown in the 
middle of Figure El When k = 4, the analogy to the smooth Schwarz P surface is 
clear. When k = 3, one can imagine a smooth embedded minimal annulus with the 
same boundary as T, and this surface is called the Schwarz H surface. 

As explained in Section ^ there are infinitely many different ways (by using 
combinations of reflections and 180° rotations that are not allowed in the smooth 
case) to extend T to a complete discrete minimal surface. Two such ways are shown 
in the upper right of Figure 21 and another two ways are shown in the center and 
right-hand side of Figure E| The two examples in Figure and the central one in 
Figure El can be extended to complete triply-periodic discrete minimal surfaces by 
180° rotations about boundary edges. The right-most example in Figure El can be 
extended to a complete triply-periodic discrete minimal surface by using horizontal 
translations perpendicular to the axis I that generate a 2-dimensional hexagonal 
grid, and then by applying vertical translations parallel to £ of length 2S(ji — jo). 
The upper-right examples in both Figures and El are applications of Method 2. 

Remark 5.1. When k — 4 and jo = —ji = 1, and when r and 5 are chosen 
properly, this T can produce the same surface as T\ produced in Subsection \5.1\ 
The way of triangulating the planar isosceles trapezoids was different in Subsection 
\5.1l but by Remark \2.1\ this is irrelevant to the minimality of the surfaces, and the 
two examples are the same in the sense that they have the same traces in K 3 . 

6. Other examples 

6.1. Discrete Schwarz CLP surface. The fundamental piece T here has six 
boundary vertices 

pi = (ar, 0, 0), P2 = (0, 0, 0), p 3 = (0, y, 0), 

Pi = (0,y,l), p 5 = (0,0,1), p 6 = (x,0,l) 
for any given fixed x, y > 0, and has one interior vertex 

p 7 = (a,b, \) . 
There are six triangles in T, which are 

(Pj,Pj+i>Pr)> j — 1 > • • • , 5 , (j>Q,pi,p 7 ) . 
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FIGURE 7. Discrete and smooth Schwarz CLP surfaces. 



The complete triply-periodic surface is generated by 180° rotations about boundary 
edges, continuing to make such rotations until the surface is complete. Every vertex 
in dT, and every vertex that is an image of a vertex in dT under these rotations, 
satisfies Equation J3J, because of the symmetry of the surface. The condition for 
Equation (J3J to hold at the interior vertex pi and all images of pj under these 
rotations is that 

2ya a a — x 

1 1 = 

^Ja 2 + \ Va 2 + (y - b) 2 ^/b 2 + {x- a) 2 

2xb b b-y = Q 

^Jb 2 + \ ^b 2 + (x- a) 2 y/a 2 + (y- b) 2 

When x = y = one explicit solution is a = b = v ^~ 1 . The fundamental piece 
T and a larger part of the resulting complete surface are shown in the left-most 
column of Figure[7|for these values of x, y, a, and b. 

For general choices of x and y, there is always a solution to the above system of 
two equations with respect to the two variables a and 6, thus giving V P7 areaT = 
0. Thus, for general x and y, the minimality condition for this example is two- 
dimensional. Fundamental pieces T for other choices of x and y are shown in the 



18 



WAYNE ROSSMAN 




FIGURE 8. Variants of the discrete Schwarz CLP surface. 



upper-center and upper-right of Figure \7\ (x = y in the center and i ^ y on the 
right). 

We can also apply Method 2 here. For example, suppose we include the reflec- 
tion of T across the plane P containing the three points (x, 0, 0), (x, 1, 0), (x, 0, 1) 
along with T to get a discrete minimal surface T\ with twelve triangles, see the 
left-hand side of Figure|Hl (Such a reflection across P would not be allowed for the 
smooth Schwarz CLP surface.) We can then extend T\ to a complete triply-periodic 
discrete minimal surface by 180° rotations about boundary edges, and this surface 
is yet another discrete superman surface. 

For a second example of applying Method 2, suppose we include the reflection 
of T\ across the plane Q containing the three points (0, y, 0), (1, y, 0), (0, y, 1) along 
with Ti to get a discrete minimal surface T2 with twenty-four triangles, see the 
right-hand side of Figure |S| (Such a reflection again would not be allowed in the 
smooth case.) We can then extend 72 to a complete triply-periodic discrete minimal 
surface by 180° rotations about boundary edges, and this gives yet another discrete 
Schwarz P surface. 

6.2. Discrete I-Wp and F-Rd surfaces. The fundamental piece T of this I-Wp 
example has six vertices 

n = (6,0,6), pa = (6,0,0), p 3 = (6,6,0), 

P4=(l, l,o), p5 = (l,a,l), Pe = prpl n WPE ■ 
There are five triangles in T, which are 

(Pl,P2,P3), (pi,P3,Pe), (P3,P4,Pe), {PA,P5,Pe), (P5,Pi,Pe)- 

By including the two images of T under the two reflections across the planes 
{(x, x, z)\x, z £ R} and {(x, y, x) | x,y S R}, we have a larger discrete surface 
with fifteen triangles. Reflecting this larger piece across all planes of the form 
{(x, y, k) I x, y £ R}, {(x, k, z) \ x,z S K}, {(k, y,z)\y,z S R} for all integers k, we 
arrive at a complete embedded triply-periodic surface in R 3 . See the right-hand 
side of Figure 

The minimality condition that Equation 10} holds at each vertex of the complete 
triply-periodic surface is 



(8) 



1 + a + a 2 - 36 - 2ab + 2b 2 = , 
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(9) a 2 + a(2 - 36) + 6 (36 - 3 + y/(l + a - 6) 2 + 2(1 - 6) 2 ) = . 

Thus, to make the surface minimal, we must find a and 6 satisfying Equations 
(JHJ-©, so the minimality condition is two-dimensional. Equation JSJ holds if 

a = - (2b - 1 + V-3 + 86 - 46 2 ) , 
and then Equation will hold if 6 satisfies 

(^3 - y/-3 + 86 - 46 2 ) (1 - b) = \/2"\/ 3 - 46 + 26 2 + + 86 - 46 2 . 

One can find such a real number 6 in a completely explicit form (although not in 
such simple forms like in Subsections 14.31 15.21 and 16.11) . 

One can similarly find a discrete analog, shown on the left-hand side of Figure EO 
of the smooth triply-periodic minimal F-Rd surface. With the simplicial structure 
chosen in Figure one can again explicitly solve the minimality condition, in the 
same way as we did for the I-Wp example. 

6.3. Trigonal example. The fundamental piece T of this H-T example has six 
boundary vertices 

a Via 1 V3 3s v^s 

Pl = ^2 , ~2~' P2 = ^2' ~2~' P3 = ( 2 ~ Y' - ^ - ' ^' 

fn 3s y/3s . .1 V3 . ,ct V3a ^ 

P4 = (2-y, — ,0), P5 = {^ 2",C), P6 = (2> 

for any given fixed 6 > 0, and has one interior vertex 

P7 = 1;{P2 +P5) ■ 
There are six triangles in T, which are 

(pj>Pj+UPr), 3 = 1 , ■ ■ ■ , 5 , (p 6 ,pi,p 7 ). 

Including the images of T under the 120° and 240° rotations about the axis 
{(0, 0, z) I z G M}, and also including the images of T and these two rotated copies 
of T under reflection across the plane {(x, y, 0) | x, y € M}, one has the larger piece 
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shown on the left-hand side of Figure ITU1 This larger piece has five boundary com- 
ponents, each contained in a plane, and these five planes bound a trigonal prism (a 
prism of height 26 over an equilateral triangle with edge- lengths 2y3). Including 
the images of this larger piece by reflecting across these five planes, and also by 
including all subsequent images of reflections across planes containing subsequent 
boundary components, one arrives at a triply-periodic discrete surface, which is em- 
bedded when a, s e (0, 1) and c G (0, 6). A larger portion of this complete discrete 
surface is shown in the central figure of Figure 1101 

The minimality condition involves three equations in the three variables a, c, s, 
and so is three-dimensional. We will not show the equations here, but they can be 
solved explicitly. For example, when 6=1, the following choices ensure minimality: 

2 + V2 3 , A h u 

a = s = , c = - (and 6=1). 

4 4 

In fact, these choices also ensure that all of the vertices of T lie in the same plane, 
and hence the fundamental piece T is planar and could be freely triangulated within 
its trace (see Remark |2.1|) . 

Furthermore, rather than using a portion of the complete surface within a trig- 
onal prism as a building block for the complete surface, one could have instead 
used a portion within a hexagonal prism as the building block. The figure on 
the right-hand side of Figure 1101 will produce exactly the same complete surface 
(again by reflecting across planes containing boundary components). In the case of 
smooth H-T surfaces, this same duality exists between building blocks in trigonal 
and hexagonal prisms, as noted in 

6.4. Discrete Fischer-Koch example. An interesting triply-periodic smooth 
embedded minimal surface was found recently by W. Fischer and E. Koch 
and is shown in the bottom row of Figure 1111 Here we give a discrete minimal 
analog of this surface, shown in the top row of Figure UTI 

The fundamental piece T of this example has eight boundary vertices 

Pl = (0,0,-1), p 2 = (0,0,-2), p 3 = (a,0,-2), p 4 = M,l), 

,a v3a „. ,a v3a „. 

= (0,0,1), p 6 = (0,0,2), ^ = (_,__, 2), P8 = (_,__,-i). 
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FIGURE 11. Discrete and smooth triply-periodic Fischer- Koch surfaces. 



for any given fixed a > 0, and has one interior vertex 



Pn 



(f4°> 



with < b < a. There are eight triangles in T, which are 

{Pj,Pj+i,P9), j = l,. ..,7, (pa,Pi,P9). 

The complete triply-periodic surface is generated by 180° rotations about boundary 
edges. In this example, the symmetry Equation Q holds at each vertex pi,...,ps in 
the resulting complete triply-periodic surface. However, getting this to hold at pg 
requires the proper choice of b. This minimality condition at pg is one-dimensional, 
and one can prove existence of a value b £ (0, a) solving it. 
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